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Abstract 


Chaos  has  been  shown  in  the  previous  project  (N00014-95-10443)  to  exhibit  in 
three-  dimensional  (3-D)  ray  tracing  for  long-range  acoustic  transmissions  in  the  ocean.  It 
imposes  a  limitation  on  our  ability  to  make  long-range  predictions.  We  have  previously 
found  that  the  sound  speed  fluctuations  in  the  upper  ocean  might  be  a  reason  for  large 
chaoticity.  Therefore,  we  inferred  that  a  smoothed  sound  speed  field  could  reduce 
chaoticity  and  thus  improve  predictability.  This  hypothesis  is  tested  and  confirmed  in  this 
project.  We  still  study  the  Heard-to- Ascension  sound  propagation,  which  has  shown  large 
chaoticity  and  extremely  limited  predictability,  in  our  previous  work,  when  NODC 
measured  sound  speed  data  were  used.  In  this  work,  the  same  NODC  data  were  smoothed 
by  using  a  sound  speed  model,  so  that  the  sound  speed  fluctuations  in  both  longitude 
direction  and  depth  direction  are  ignored.  We  use  the  smoothed  sound  speed  to  perform 
3-D  ray  tracing.  Numerical  results  show  that  chaoticity  was  significantly  reduced  and 
consequently  predictability  greatly  improved.  The  ray  paths  show  a  pattern  similar  to 
those  constructed  using  the  original  NODC  data.  The  travel  times,  however,  were  less 
than  the  previous  results  by  about  10  seconds,  with  a  relative  error  of  0.  16%.  This  error 
was  produced  by  the  smoothed  sound  speed. 
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1.  Introduction 


Chaos  is  defined  to  be  aperiodic  bounded  dynamics  in  a  deterministic  system  with 
sensitive  dependence  on  initial  conditions.1  In  1988,  Palmer  et  al.2  demonstrated  that 
chaos  appeared  in  two-dimensional  (2-D)  acoustic  ray  tracing  in  the  ocean.  It  is  then 
called  “ray  chaos”  by  the  ocean  acoustics  community  and  further  studies  3‘4’5’6’7’8’9’10’11 
have  been  performed  following  the  Palmer  et  al.’s  work. 

These  previous  studies  are  important.  However,  they  all  deal  with  2-D  ray 
tracing.  The  2-D  equations  are  not  valid  for  long-range  transmissions,  while  ray  chaos  is 
a  phenomenon  in  long-range  sound  propagation.  In  other  words,  the  2-D  equations  may 
lose  validity  before  chaos  is  manifest.  This  can  be  seen  from  Smith  et  al.’s  work.8  In  that 
work,  a  predictability  horizon  of  1000  to  2000  km  has  been  estimated  by  using  Lyapunov 
exponents,  but  the  ray  equations  used  in  that  work  are  only  valid  for  about  100  km.  This 
implies  that  the  predictability  was  limited  by  the  validity  of  the  2-D  ray  equations,  rather 
than  chaos. 

For  long-range  transmissions,  we  currently  have  three  sets  of  ray  equations.  They 
are  horizontal  ray  equations,  ’  Hamiltonian  equations  in  HARPO,  and  3-D  ray 
equations  in  ellipsoidal  coordinates.15  These  equations  are  valid  for  very  long  ranges. 
Chaos  may  be  a  major  factor  that  determines  the  predictability  of  these  ray  models. 
Therefore,  we  need  to  used  these  models  to  study  chaos  to  understand  our  ability  to  make 
long-range  predictions.  We  have  recently  studied  chaos  in  horizontal  ray  tracing, 16and 
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then  in  the  previous  project  (ONR  Grant  N000 14-95- 10443),  we  have  studied  chaos  in  3- 
D  ray  tracing.17’18’19 

We  have  studied  two  3-D  cases  in  the  previous  work.17,19  The  first  is  the  Heard-to- 
Ascension  (H-A)  sound  propagation,  and  the  second  is  the  Califomia-to-Hawaii  (C-H) 
propagation.  Chaos  appeared  in  both  cases.  However,  there  was  a  significant  difference 
between  the  two  cases.  In  the  H-A  case,  chaoticity  was  very  strong,  and  predictability 
was  extremely  limited  by  chaos,  In  the  C-H  case,  however,  chaos  was  weak,  and 
predictions  were  made  with  high  precision.  Comparison  between  the  two  cases  suggests 
that  there  is  an  apparent  difference  in  sound  speed  data,  besides  the  difference  in 
propagation  range.  The  sound  speed  data  used  for  the  H-A  sound  propagation  have  large 
fluctuation  in  the  upper  ocean,  as  can  be  seen  from  Fig.  1.  The  sound  speed  fluctuation  in 
the  C-H  case,  however,  is  very  small,  as  seen  from  Fig.  2.  Therefore,  the  smoothness  of 
the  sound  speed  data  in  the  C-H  case  might  be  a  reason  for  the  weak  chaos.  This  leads  to 
a  hypothesis:  a  smoothed  sound  speed  field  can  reduce  chaoticity  and  thus  enhance  our 
ability  to  make  long-range  predictions.17 19 

The  objective  of  this  project  is  to  test  this  hypothesis.  To  study  the  feasibility  of 
using  smoothed  sound  speed  data  to  reduce  chaoticity,  we  tackle  the  H-A  problem  that 
had  very  large  chaoticity  in  our  previous  work.1719  We  use  the  same  sound  speed  data  as 
those  used  in  our  previous  work,  but  before  ray-tracing  we  fit  the  data  to  a  mathematical 
model  and  use  the  model  to  provide  a  smoothed  sound  speed  field.  This  is  described  in 
section  two.  In  section  three,  we  use  the  smoothed  sound  speed  data  to  perform  ray 
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tracing,  and  use  Lyapunov  exponent  to  quantify  chaoticity.  The  numerical  results  will  be 
compared  with  those  of  our  previous  work  to  test  the  feasibility  of  the  proposed  method. 
In  section  four,  presented  are  summary  and  discussions.  We  focus  on  the  correctness  of 
the  ray  tracing,  considering  that  the  smoothed  sound  speed  field  might  lose  reality,  which 
could  in  turn  affect  the  correctness  of  ray  tracing.  We  wonder  if  there  is  a  compromise 
between  an  improved  predictability  and  the  correctness  of  ray-tracing  results. 

2.  Smoothing  Sound  Speed  Data  with  Munk’s  Model 

2. 1  Characteristics  of  Sound  Speed  Profiles 

From  the  sound  speed  data  provided  by  the  National  Oceanic  Data  Center 
(NODC),  we  can  see  that 

•  the  sound  speed  fluctuates  at  the  upper  ocean,  and 

•  for  a  given  latitude,  the  sound  speed  profiles  (SSP)  at  different  longitude  are 
similar  in  shape. 

These  characteristics  can  be  seen,  for  example,  from  Fig.  1 .  We  have  studied  all  the 
NODC  data  for  the  Southern  Atlantic  Ocean,  and  this  is  true  for  all  the  data  files  we  have 
received  from  NODC,  as  depicted  in  the  Appendix. 

The  first  characteristic  might  be  a  major  factor  that  causes  chaos  in  ray  tracing. 
The  second  characteristic,  however,  may  lead  to  a  method  to  reduce  chaoticity.  It  allows 
us  to  use  a  model  to  approximate  all  the  SSP  for  a  given  latitude.  In  this  way,  we  can 
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smooth  the  NODC  sound  speed  data  to  reduce  chaoticity.  Such  a  model  is  implemented 
as  follows. 

2.2  Fitting  Data  to  Munk’s  Canonical  Model 

In  ocean  acoustics,  Munk20  has  presented  a  sound  speed  model  expressed  as 

c(r)  =  c«[l  +  £(77  +  e~,?  “l)]>  (*) 

where 

i)  =  l(r  -  ra)  /  B.  (2) 

Here  c  is  the  sound  speed;  r  is  the  depth;  ca  and  ra  are  the  sound  speed  and  the  depth  at 
channel  axis,  respectively;  B  is  the  scale  depth;  sis  the  perturbation  coefficient;  and  77  is  a 
dimensionless  distance.  This  model  is  commonly  considered  as  a  model  for  deep  ocean. 
By  plotting  the  NODC  data,  however,  we  note  that  for  a  given  latitude,  all  the  sound 
speed  profiles  including  shallow  water  SSP  are  similar  in  shape,  except  that  they 
terminate  at  different  depths.  This  suggests  that  we  can  use  Munk’s  model  to 
approximate  all  the  SSP  for  a  given  latitude.  This  can  be  done  by  fitting  the  NODC 
measured  sound  speed  data  to  Munk’s  model,  described  as  follows. 

Let  c  be  the  sound  speed  computed  using  Munk’s  model,  and  C  be  the  NODC 
measured  data.  Then,  the  difference  between  the  model  output  and  die  measured  data  can 
be  quantified  using 

(3) 

<=1  7=1 
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Equation  (3)  means  that  we  have  N  SSP  for  a  given  latitude,  and  each  SSP  is  sampled  at  n 
depth  points.  Here,  the  depth  is  not  equally  discretized.  The  NODC  data  have  the 
standard  depths  listed  in  Table  1 .  In  Eq.  (3),  x  represents  the  model  parameters  in  Eq.  (1) 
and  (2),  i.e.,  ca ,  s,  i),  ra ,  and  B. 

Now  we  can  fit  the  NODC  measured  data  to  Munk’s  model  by  solving  the 
following  problem: 

Minimize  {M(x)},  (4) 

subject  to  the  constraints 

c1  <  c  <cu . 

el  <s<su, 

rjl  <rt<7f ,  (5) 

rla  <ra  <r“ , 

Bl  <B<BU, 

where  superscript  /  represents  the  lower  bound  of  the  parameter,  and  u  the  upper  bound. 
In  the  following  calculations,  we  use  0  <  s  <  0.01,  and  0.4  <  B  <  1.7 .  The  ranges  of 
ca  and  ra  vary  substantially  with  latitudes,  and  thus  are  determined  directly  from  the  plots 
of  the  measured  NODC  data  in  the  Appendix. 

Problem  (4)  is  a  constrained  optimization  problem,  and  can  be  solved  using  the 
Complex  method.21 ,22  Using  this  method,  we  can  estimate  the  model  parameters  in 
Eqs.(l)  and  (2).  We  have  NODC  measured  sound  speed  data  of  the  Southern  Atlantic 
Ocean.  The  latitude  of  the  data  ranges  from  8°S  to  52°S  with  1°  increment.  We  first  sort 
the  NODC  data  into  numerical  orders  of  latitudes,  and  get  45  data  sets  corresponding  to 
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45  latitudes.  Then  we  fit  each  data  set  in  Munk’s  model  (1)  using  the  Complex  method 
of  optimization.  The  model  parameters  estimated  are  listed  in  Table  2. 

Figure  3  shows  typical  results  of  fitting  the  NODC  data  to  Munk’s  model.  We 
can  see  that  the  model  predictions  are  in  good  agreement  with  the  NODC  measured  data, 
not  only  at  low  latitude  (Fig.  3-1),  but  also  at  high  latitude  (Fig.  3-2).  For  mid-latitudes, 
however,  the  discrepancy  between  the  NODC  data  and  the  model  predictions  is  larger 
(Fig.  3-3).  This  is  because  that  the  measured  data  fluctuated  substantially  at  these 
latitudes.  The  output  of  Munk’s  model  is  still  a  good  approximation  as  an  averaged  SSP. 

2.3  A  Sound  Speed  Model  of  Southern  Atlantic  Ocean 

By  fitting  the  NODC  measured  sound  speed  data  to  Munk’s  model,  we  can  now 
use  Eq.  (1)  to  approximate  all  the  SSP  for  a  given  latitude.  In  this  way,  we  can  smooth 
the  sound  speed  in  both  the  longitude  direction  and  the  depth  direction.  Suppose  that  the 
sound  speed  data  are  sampled  at  N  different  latitudes,  expressed  as  </\ ,  . . . ,  $N .  Then  we 

have  N sound  speed  profiles,  c(r)i,c(r)2,...,c(r)N ,  which  can  be  approximated  using  Eq. 
(1).  For  a  latitude  located  between  <j>t  and  <f>M ,  the  sound  speed  can  be  calculated  using 
the  linear  interpolation  between  c{r)i  and  c(r)M ,  i.e., 

=  c(r),  + — - - - {.<!>- (pt)  (6) 

Yi+\  Yi 

Equation  (6)  is  a  sound  speed  model  we  shall  use  to  perform  ray  tracing  for  Heard-to- 
Ascension  sound  propagation.  This  model  ignores  the  sound  speed  fluctuations  in  both 
the  longitude  and  the  depth  direction.  It  might  be  a  good  approximation  considering  that 
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these  fluctuations  are  relatively  small  compared  with  that  in  the  latitude  direction.  We 
shall  investigate  the  feasibility  of  using  this  model  to  improve  the  predictability  of  ray 
tracing,  and  discuss  the  correctness  of  this  model  in  section  four. 


3.  Numerical  Ray  Tracing 

We  have  just  formulated  a  sound  speed  model  for  the  Southern  Atlantic  Ocean. 
The  model  smoothes  the  NODC  measured  sound  speed  data  in  both  the  longitude 
direction  and  the  depth  direction.  Now  we  use  this  model  for  ray  tracing.  We  have  two 
objectives.  First,  we  test  the  hypothesis:  a  smoothed  sound  speed  field  can  reduce 
chaoticity,  and  improve  predictability.  Second,  we  examine  the  correctness  of  the  ray 
tracing — we  wonder  if  the  smoothed  sound  speed  field  is  a  good  approximation  to  reality. 


3.1  Methods 


For  comparison,  we  still  use  the  3-D  ray  equations  that  we  used  in  our  previous 
work  (17,19).  They  are 


d(f>  _  cos  9  cos  a 
ds  ju-r 

dX  cos#sina 


ds  (v-r)cos^’ 
dr 


ds 


=  sin#, 


da  costftandsina  .  _  .  (  1 

—  = - - - +  sin  6/ sin  a  cos  a 

ds 


v-r 
sin  ad 


1 


A 


cos  a 


\v-r  fj.-rj 
In N 


H~r  dtj>  (v-r)cos^  dX)  cos# 


(7) 

(8) 

(9) 


(10) 


and 
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to  estimate  propagation  range.  The  derivation  of  equation  (15)  has  previously  been 
presented  in  Ref.(17,19). 
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These  equations  are  numerically  integrated  by  using  die  fourth-order  Runger- 
Kutta  method  with  adaptive  step  size  control,24  the  same  as  that  in  our  previous  work. 

The  calculation  accuracy  is  controlled  by  specifying  a  maximum  fractional  error  (10"6 )  in 
any  single  integration  step.  Sound  speed  is  estimated  using  the  model  described  in  last 
section.  An  eigenray  is  constructed  in  the  following  way.  First,  we  specify  a  launch 
grazing  angle  (0t)  and  search  for  the  launch  azimuth  (  a, )  that  drives  the  ray  to  pass 
within  a  given  error  tolerance  for  both  the  latitude  and  the  longitude  of  the  receiver. 
Successive  Shooting  method25  is  used  for  the  searching  process.  Then,  we  change  0,  and 
again  search  for  a  corresponding  at .  The  searching  process  terminates  if  a  ray  passes 
within  the  given  tolerances  for  the  receiver’s  coordinates  (</>r,Ar,rr ).  The  locations  of  the 
source  and  the  receiver  in  this  calculation  are  listed  in  Table  3.  These  data  are  obtained 
from  the  published  references  26,27  The  receiver’s  parameters  are  those  of  hydrophone  23 
located  south  of  the  Ascension  Island.27 

Chaos  is  measured  by  using  Lyapunov  exponent.  Any  chaotic  system  must  have 
at  least  one  positive  Lyapunov  exponent,  and  the  magnitude  of  the  exponent  is  a  measure 
of  chaoticity.  In  the  following,  we  follow  Wolf  et  al.’s  definition.  The  z'th  one- 

9ft 

dimensional  Lyapunov  exponent  is  defined  as 

4  =  lim  ^  log2  ^4^ ,  (16) 

t  a(°) 

where  p  is  the  length  of  the  ellipsoidal  principle  axis,  and  h  are  ordered  from  the  largest 
to  the  smallest.  Wolf  et  al.’s  method  is  used  to  estimate  the  Lyapunov  exponent,  and 
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their  published  FORTRAN  code28  is  adapted  for  the  calculations  in  this  work.  Along  a 
ray  path,  we  numerically  integrate  the  variational  equations  given  in  Ref.  (17),  and  then 
use  Wolf  et  al.’s  algorithm  to  estimate p  in  Eq.  (16)  and  compute  Lyapunov  exponents. 
The  magnitude  of  the  largest  Lyapunov  exponent  of  an  eigenray  is  used  to  quantify  chaos. 

3.2  Chaoticity 

Figure  4  shows  the  vertical  ray  paths  of  eigenrays  constructed  using  the  smoothed 
sound  speed.  The  corresponding  ray  parameters  are  listed  in  Table  4.  Now  we  quantify 
chaoticity  using  the  magnitude  of  the  largest  Lyapunov  exponent  of  an  eigenray.  Using 
the  method  described  in  section  3.1,  we  estimate  the  Lyapunov  exponents  for  all  the 
eigenrays  in  Fig.  4.  The  largest  Lyapunov  exponent  for  each  eigenray  is  given  in  Table  5. 
For  comparison,  Table  6  gives  the  largest  Lyapunov  exponent  for  the  eigenrays  computed 
using  the  original  NODC  data  in  our  previous  work  (17,19). 

Comparing  Table  5  with  Table  6,  we  can  see  that  the  eigenrays  constructed  using 
the  smoothed  sound  speed  have  apparently  smaller  Lyapunov  exponents  than  those 
computed  using  the  original  NODC  sound  speed  data.  This  suggests  that  using  the 
smoothed  sound  speed  significantly  reduced  chaoticity.  This  is  consistent  with  our 
hypothesis. 

We  can  see  from  Table  6  that  the  near-axial  ray  paths  (indicated  by  small 
identification  number)  had  larger  Lyapunov  exponents.  This  is  because  that  these  rays 
were  constructed  using  the  original  NODC  data  that  have  large  fluctuation  near  the 
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channel  axis  (see  Fig.  1).  In  Table  5,  however,  the  near-axial  rays  have  smaller  Lyapunov 
exponents,  because  these  rays  are  computed  using  the  smoothed  sound  speed  that  ignores 
the  large  fluctuation  in  the  upper  ocean.  Comparison  between  Table  5  and  Table  6 
suggests  that  the  large  fluctuation  in  the  upper  ocean  is  a  main  cause  for  the  chaos  in  ray 
tracing  and  the  smoothed  sound  speed  reduces  chaoticity  by  ignoring  this  fluctuation. 

3.3  Predictability 

We  have  just  shown  that  the  chaoticity  was  significantly  reduced  by  using  the 
smoothed  sound  speed.  As  a  result,  the  predictability  will  be  significantly  improved,  as 
shown  in  this  section. 

Recall  that  in  our  previous  work  (17,19),  we  used  original  NODC  sound  speed 
data  to  construct  the  eigenrays  between  Heard  Island  and  Ascension  Island.  The 
prediction  accuracy  was  typically  less  than  0.005  degree  in  the  receiver’s  longitude,  and 
less  than  1  km  in  the  receiver’s  depth. 

Compared  with  this  result,  the  prediction  accuracy  of  this  work  was  significantly 
improved.  This  can  be  seen  from  Table  4.  From  this  Table,  we  can  see  that  the  first  four 
eigenrays  had  the  error  of  less  than  0.000005  degree  in  the  receiver’s  longitude,  and  less 
than  0.0005  km  in  the  receiver’s  depth.  Although  eigenray  5  through  7  had  relatively 
large  errors.  The  prediction  accuracy  of  these  rays  was  still  better  than  that  of  our 
previous  result. 
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4.  Summary  and  Discussions 

The  objective  of  this  report  is  to  investigate  the  feasibility  of  using  smoothed 
sound  speed  to  reduce  chaoticity  in  acoustic  ray  tracing  in  the  ocean.  This  is  done  by  (1) 
fitting  the  NODC  measured  sound  speed  data  to  Munk’s  model,  (2)  performing  ray 
tracing  with  the  model-produced  sound  speed,  and  (3)  comparing  chaoticity  with  that 
previously  produced  using  original  NODC  data.  Numerical  results  show  that 

1 .  Chaoticity  was  significantly  reduced  by  using  the  smoothed  sound  speed  and 
consequently,  the  predictability  was  greatly  improved. 

2.  Large  fluctuation  of  sound  speed  in  the  upper  ocean  is  a  main  cause  for  the  chaos 
in  acoustic  ray  tracing  in  the  ocean. 

Now  we  have  a  question:  Did  we  lose  the  correctness  of  the  ray  tracing  while 
enjoying  the  greatly  improved  predictability?  Note  that  we  used  a  model  to  smooth  the 
NODC  data.  The  model’s  output  is  just  an  averaged  SSP  for  a  given  latitude.  This  is  an 
approximation.  There  must  be  some  discrepancy  between  the  result  calculated  using  the 
original  NODC  data  and  the  smoothed  data.  To  examine  the  correctness  of  the  ray 
tracing  with  the  smoothed  sound  speed,  let  us  compare  the  results  with  those  computed 
using  the  original  NODC  data. 
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First,  let  us  compare  the  ray  paths.  According  to  our  previous  results  (17,19),  the 
vertical  paths  started  in  a  shallow  sound  channel  near  Heard  Island.  This  pattern  persists 
along  a  distance  of  about  3000  km.  Then  a  transition  began  from  the  shallow  channel  to  a 
deep  channel.  The  rays  entered  the  deep  channel  at  the  propagation  range  of  about  4000 
km,  and  finally  reached  Ascension  Island.  This  is  consistent  with  the  ray  paths 
constructed  using  the  smoothed  sound  speed  (see  Fig.  4). 

Then,  let  us  compare  the  travel  times.  The  largest  travel  time  computed  using  the 
smoothed  sound  speed  was  about  6248  s  (see  Table  4),  while  that  produced  using  the 
original  NODC  sound  speed  data  was  about  6258  s.  There  was  10-second  discrepancy 
between  the  two  results.  In  other  words,  using  the  smoothed  sound  speed  to  approximate 
the  original  NODC  data  led  to  a  relative  error  of  0.16%  (10/6258)  in  ray  tracing.  This 
small  error  might  be  acceptable  for  such  a  long-range  (9223  km)  sound  propagation.  On 
the  other  hand,  it  was  impossible  to  get  an  accurate  3-D  eigenray  using  the  original 
NODC  data  as  mentioned  in  our  previous  work.  Using  the  smoothed  sound  speed  made 
it  possible,  however. 
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Table  1  NODC  standard  depth 
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Table  2  Fitting  NODC  data  to  Munk’s  Model 


Model  Parameters 


Latitude 

(degree) 


8S 

9S 

10S 

11S 

12S 

13S 

14S 

15S 

16S 

17S 

18S 

19S 

20S 

21S 

22S 

23  S 

24S 

25  S 

26S 

27S 

28S 

29S 

30S 

31S 

32S 

33  S 

34S 

35S 

36S 

37S 

38S 

39S 

40S 

41S 

42S 

43  S 

44S 

45S 

46S 


Ca 


1.483221 
1.482676 
1.482837 
1.482376 
1.483215 
1.482321 
1.481682 
1.481949 
1.482206 
1.482156 
1.481896 
1.482056 
1.482463 
1.482562 
1.482844 
1.482559 
1.482920 
1.483193 
1.481971 
1.482843 
1.482556 
1.48165 6 
1.482742 
1.484858 
1.484963 
1.48499 9 
1.484671 
1.484995 
1.484234 
1.482627 
1.488353 
1.487902 
1.485842 
1.481831 
1.482055 
1.480181 
1.477453 
1.474089 
1.471392 


0.002194 

0.002421 

0.002226 

0.002211 

0.002044 

0.002610 

0.003434 

0.003451 

0.003358 

0.003563 

0.004038 

0.004727 

0.004605 

0.005757 

0.005775 

0.005437 

0.005638 

0.006417 

0.006695 

0.007946 

0.008726 

0.006837 

0.007310 

0.009990 

0.009949 

0.009894 

0.009974 

0.007662 

0.009995 

0.007826 

0.009999 

0.009176 

0.005718 

0.004701 

0.009847 

0.004828 

0.093715 

0.003699 

0.009948 


0.732704 

0.747110 

0.715887 

0.715331 

0.700050 

0.747294 

0.818365 

0.835788 

0.812779 

0.848575 

0.896096 

0.913893 

0.927040 

0.997749 

0.996901 

0.999972 

0.999393 

0.999880 

0.983354 

1.057264 

1.090114 

1.012344 

1.097952 

1.176254 

1.199543 

1.195681 

1.132879 

1.099913 

1.191613 

1.080397 

1.278593 

1.199978 

1.038763 

0.883005 

0.991699 

0.796828 

0.689430 

0.484000 

0.380606 


B 


0.472351 

0.498450 

0.464047 

0.456252 

0.452106 

0.511343 

0.612715 

0.619859 

0.602596 

0.645464 

0.706073 

0.782972 

0.780659 

0.923985 

0.924338 

0.890491 

0.917236 

0.987625 

1.020320 

1.159711 

1.240750 

1.029022 

1.136131 

1.414319 

1.451312 

1.448759 

1.459401 

1.218559 

1.431734 

1.181853 

1.634727 

1.481553 

1.081195 

0.870412 

1.469271 

0.896648 

1.315256 

0.702803 

1.384595 


Standard 

Deviation 


0.00829 

0.00559 

0.00346 

0.00460 

0.00297 

0.00379 

0.00439 

0.00369 

0.00064 

0.00065 

0.00023 

0.00389 

0.00429 

0.00442 

0.00280 

0.00384 

0.00549 

0.00141 

0.00947 

0.00225 

0.00907 

0.00290 

0.00382 

0.00656 

0.00421 

0.00553 

0.00095 

0.01729 

0.01088 

0.00775 

0.03697 

0.03103 

0.02093 

0.00542 

0.02515 

0.01789 

0.00232 

0.00126 

0.00509 
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47S 

1.467830 

48S 

1.464943 

49S 

1.462756 

50S 

1.460705 

51S 

1.457501 

52S 

1.460736 

0.002195  0.300742 

0.002192  0.066605 

0.002269  0.143387 

0.002239  0.010283 

0.002644  0.001047 

0.002236  0.052927 


0.400003  0.00440 

0.419463  0.00163 

0.401664  0.00208 

0.400460  0.00241 

0.401966  0.00031 

0.400119  0.00327 
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Table  3  Parameters  of  source  and  receiver 


Source 

Receiver 

Latitude  (Deg) 

53°22’S 

8°4.2’S 

Longitude  (Deg) 

74°30’E 

14°25.2’W 

Depth  (km) 

0.175 

0.832 

Table  4  Eigenray  parameters 


No. 

Launch  angles  (degree) 

Travel 

time(s) 

Range 

(km) 

Er* 

Ex** 

Grazing 

Azimuth 

1 

- 1 .425323486328 1 25  265 .6 1 34263277054 

6248.5 

9222.4 

0.0004 

3  x  10"6 

2 

-1.639988307326689  265.6392212957143 

6248.5 

9222.4 

0.0005 

6  x  10~7 

3 

2.2421875 

265.724887589873 

6248.6 

9222.4 

0.0002 

lxKT6 

4 

4.091796875 

265.903826713562 

6248.1 

9222.3 

0.0004 

4  x  10-6 

5 

-4.2421875 

265.9238282442446 

6248.1 

9222.3 

0.001 

7  x  10-4 

6 

-4.60302734375 

265.9766235686838 

6247.7 

9222.2 

0.0003 

7 

0 

i-H 

X 

7 

6.467772959697868 

266.243524546997 

6245.6 

9222.2 

0.0001 

2  x  10-4 

*  Er  =  \rend  -  0.832| 

**E,  =  \^end  “  0.000005) . 
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Table  5  Lyapunov  exponent  of  eigenrays 


Eigenray  No. 

Largest  Lyapunov  exponent 
(1/km) 

1 

2.6  xlO"3 

2 

2.8  x  10“3 

3 

4.2  x  10'3 

4 

1.7  xlO-2 

5 

1.7  xlO-2 

6 

1.8  xlO"2 

7 

2.0  xlO"2 
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Table  6  Lyapunov  exponent  of  eigenrays  computed  using  NODC  data 


Eigenray  No. 

Largest  Lyapunov  exponent 
(1/km) 

1 

1.1  x  lO'1 

2 

1.1  xlO-1 

3 

3.8  xlO"2 

4 

4.0  xlO”2 
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Depth  (m) 


Fig.  1  NODC  sound  speed  profiles  of  the  Atlantic  Ocean  (Latitude  15  °S,  Longitude  from 
1W  to  13  CW).  The  sound  speed fluctuates  significantly  at  the  upper  ocean,  but  they  are 
similar  in  shape 


Fig.  2  Levitus  sound  speed  profiles  of  Pacific  Ocean  (Latitude  25.5  °N,  Longitude  from 
159.5  °W  to  120.5  °W).  These  profiles  are  relatively  smooth  compared  with  the  NODC 
data  in  Fig.  1 
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3-3  Latitude:  36S 


Fig.  3  Typical  results  of  fitting  NODC  data  to  Munk’s  model 
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Fig.  4  Vertical  ray  paths  of  eigenrays  between  Heard  Island  and  Ascension  Island 
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Appendix 

Fitting  NODC  sound  speed  data  to  Munk’s  model 

In  the  following  figures,  the  numbers  in  the  legend  denote  longitudes,  and  the  model 

output  is  a  thick  line. 
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